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Abstract 

Limit tlieorems are proved for quadratic forms of Gaussian random 
fields in presence of long memory. We obtain a non central limit theo- 
rem under a minimal integrability condition, which allows isot ropic and 
anisot ropic models. We apply our limit theorems and those of iGinovianI 
l|l999l l to obtain the asymptotic behavior of the empirical covariances of 
Gaussian fields, which is a particular example of quadratic forms. We 
show that it is possible to obtain a Gaussian limit when the spectral den- 
sity is not in . Therefore the dichotomy observed in dimension d = 1 
between central and non central limit theorems cannot be stated so easily 
due to possible anisotropic strong dependence in d > 1. 
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1 Introduction 

Many statistical methods proposed for long memory proces ses are based on 
statistics, which can be expressed using quadratic fo rms (see |Bera4 Il994 for 
a review). The empirical covariance sequence (see iHoskine . 



199(: ; Hannaril . 



19761) ■ the estimation of the long memory parameter using Whittle con trast (see 



Fox and Taggu . T986I : Giraitis and Surgailid . ll990l : Giraitis and Taggu . 1999L for 
exam ple) or using the integrated periodogram (see lLobato I. and Robinson P.M.I . 
199i), some change point detection procedure (see lBeran and Terrinl . ll996[ ) rely 
on quadratic forms. 

It would be interesting to validate the same statistical tools for the study of 
random fields having long memory. As a first step, we study in this paper the 
convergence of the quadratic forms in the d— dimensional (d > 1) case. 

Let us introduce our framework. Let X — (X„)„gg£i be a stationary 
random field having long memory, i.e. its covariance sequence {r{n))ni=z'^ is not 
a summable series. We assume that X is a Gaussian process and admits the 
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linear representation 



Xn = / a(a;)e* 

J E 



(1) 



where W is the Gaussian white noise spectral measure and where the function 
a is in L?'{E). < •, • > denotes the inner product in W^. Hereafter, we denote 
by / the spectral density of X, f is proportional to |a(a::)p, and we have 



rin) = fn 



„'i<n,A> 



/(A)dA. 



The quadratic forms associated to X and {gi)i^fid e ^ (N ) are defined as 



gi-jXiXj 



(2) 



(3) 



ieA„ jeA„ 



where An = {1, . . . , n}'^. The statistics Qn can be rewritten as a functional of 
the periodogram /„ of X. Indeed let 



/„(t) = (2™)-'* >: XkXie'<''-'^'> 



we have 



(27rf / g{t)In{t)dt, 



(4) 



-i<j,t> 



where 

g{t) = (27r)-'^ ^ 

In dimension one (d = 1), many papers deal with the asymptotic behav- 
ior of Q„ for long memory processes. When the intensity of the memory is 
not too strong or when specific conditions are imposed to g in order to kill 
the effect of the long memory involved by /, the convergence rate and the 
asymptotic normality obtained in short memory are pre s erved . The asymptotic 
norma lity of the quadratic forms i s prov ed by lAvraml ( 19881 ). I Fox and Tag qui 
(IiQBTI) and Ginovian and Sahakian ( 2005 ) in the Gaussian case and extended by 
iGiraitis and Surg aihd l|l99(lh to include the non-Gaussian linear case. Avram's 
result is restricte d to / and q in L^, bu t requires only integrabil i ty con ditions. 
On the opposite, Fox and TaqqrJ ( 1987 ). Ginovian and SahakianI ( 20051 ) specify 
th e behavior of / and g a round zero. The extension to linear m odels obtained 
by Giraitis and Surgailij 1 1990f ) holds under the assumptions of Avram (1988) 
or lFox and TaqquI (|l987[ ). 

Under o ther conditions on / and g fixing the b e havior of both functions 
around zero, Rosenblatt ( 1961 ). Fox and Taqqu (1985), Terr in and Taggu ( 1990l ) 
prove non central limit theorems for Gaussian processes. The normalization is 
not standard (i.e. it is not equal to n^'^/^) and the Umiting process is no t Gaus- 
sian. The extension to linear models is proved in IGiraitis et al (Il998h in the 
general context of bivariate Appell polynomials. 
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Fixing the behavior of / and g around zero, the dichotomy in d = 1 can be 
summarize as follows. Assume that 



f{x) ^ |a;p" when 2: — >■ 0, 
g{x) ^ \x\^^ when a; — s> 0. 

The limiting distribution depends on the value of (a + /?). 

• Ifa + /3> — then the limit is Gaussian and the normalization is n~^/^ 
{Qn satisfies a central limit theorem). 

• li a-\- B < — then the limiting distribution is not Gaussian and the nor- 

4 ° 

malization in this case is which is lower than n~'^/'^ {Qn satisfies 

a non central limit theorem). 



The case a + /? = - is investigated in Ginovian and SahakianI ( 20051 ) (up to a 

slowly varying function in / and g which vanishes at zero). In this case, the 
limit remains gaussian. 

In the d-dimens ional case (d > 1), a ce ntral lim it theorem has been proved 



3 ( I999I ). in the same spirit as lAvraiii (fl988i) (i.e. / and g are in 



by iGinoviar 

and satisfy an integrability condition). This result is recalled in Section [3TTJ 
When f or g are not in LP' , some central limit theorems have be en proved under 
some specific conditions o n these funct i ons, s ee Theorem 3.1 in Doukhan et al.l 
( 19961 ) and Lemma 3.3 in lBoissv et"al] (|2005l ). 

On the other hand, few results are availab le on non-central limit theorems 
in dimension d > 1. In Doukhan et al. ( 19961 ). th e asymptot i c of Q n is inves- 
tigated when the memory of X is "isotropic" (see Lavancien ( 2005 ) Definition 
1), i.e. the memory of X is due to in the expression of / where | • | de- 

notes the Euclidian norm. The a uthors obtain a non-central limit theorem very 
similar to IFox and Taqqui (1985). However the investigation of "anisotropic" 



memory and more general functions g (including the particular case of the em- 
pirica l covariance) is of great interest for statistical applications. In iLavancier 
( 2005h . different classes of random fields having long memory are described: The 
important role of anisotropic processes in terms of modeling is illustrated. 

Our paper is organized as foll ows. In S ection [21 after a reminder of the 



central limit theorem obtained by Ginovian ((l999l ). we present a non-central 
limit theorem under a general integrability condition on / and g. We show that 
this condition is satisfied by a large class of isotropic and anisotropic models 
(including all previous results). Since random fields in d > 1 may exhibit much 
more different types of long memory than processes \i\ d — 1, the dichotomy 
between central and non central limit t heorems is not so simp l e to fo rmulate. As 
already highlighted in Remark 4.2 of Dobrushin and Maior ( 1979I ). a spectral 
singularity outside zero may lead to new limiting results. This fact is illustrated 
through examples in Section 12.31 
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In Section [21 we apply limit theorems for quadratic forms to the convergence 
of the empirical covariance series. Contrary to the classical setting in dimension 
one, we show that it is possible to obtain a Gaussian limit in dimension d > 1 
when / is not in L^. 



2 Convergence of quadratic forms 



2.1 Central limit theorem 

We recall in this section the centra l limit theorem o btained by AvramI (1988) in 
dimension d = I and extended by iGinovia to d > 1. This result will 

be useful in section |31 where the asymptotic behavior of the covariance series is 
investigated. 



Theorem 1 (jAvraml (|l988l ). [Cinovianl (|l999l )). Let X be a random field as in 
([Ip and denote f its spectral density. Let Qn, defined in (Q), be the quadratic 
form associated with g and X . If f € L^ and g G L"^ with p > 2, q > 2 and 
l/p+l/q < 1/2, then 



n^/^Qn - E{Qn)) ^ M (0,2{27rf'' / f\t)g\t)dt 



(5) 



where — > denotes the convergence in law. 

This result provides a large setting since no specific form for functions / and 
g is required, but only integrability conditions. In re t urn, / and g are restricted 
to L^. A similar result is obtained in lDoukhan et al.l (jl996l ). Th. 6.2., but under 
semi-parametric assumptions on / and g. 

When f oi g are not in L^ but a com pensation betwe e n the m occurs, a 
central limit theorem may still arise as in IFox and TaqquI (jl985l ) for d = 1. 
To our knowledge, two situations in dimension d > 1 have been investigated 
so far. In di mension d = 2, ass uming / and g have a tensorial product form 



around zero. iBoissv et al~ ( 2005 ) obtain a central limit theorem , which can be 
straightforwardly extended to d > 2. In iDoukhan etlU] (|l996[ ). Th. 6.1, the 
authors prove the same result when / has an isotropic form around zero and 
when g is regular enough and g{0) — 0. 

It is known that the lonely condition < f^it)g^(t)dt < cxj is not sufficien t 
to obtain a central limit theorem, as proved in lGinovian and" SahakianI (l2005h . 
Prop. 2.2. A general setting leading to a central limit theorem (included / or g 
not in L^) is still missing. This issue is not investigated in this paper. 



2.2 Main result : non-central limit theorem 

The following theorem gives the convergence of the quadratic forms under the 
general condition (H) below on / and g. Contrary to the previous studies in 
the one dimension case or in dimension d > I, this condition does not require to 
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specify the behavior of / around its singularities. Such a condition is satisfied 
by many examples as shown in Section [2.31 

Let ^ be a function defined on R'', hereafter we denote by Qp) (p G N) 
the periodic function with period 27rp (with respect to each component), that 
coincides with ( on the set pE"^, 



(6) 



Theorem 2. Let X be the linear field defined by (QP and Qn the quadratic form 
defined by 

Assume that for all x ^ E a(— x) = a{x), and g{—x) — g{x). Moreover, 
assume that 

\a{x) = a(i)(x)Li(x), 
\g{x) = g(i){x)L2{x), 

where a is a homogeneous function of degree a (i.e. Vc > 0, Va: G E, d{cx) = 
c"a(x)J, where g is a homogeneous function of degree 2/3 and where Li and 
are bounded functions, continuous at zero and Li(0) ^0 (i = 1,2). 
If the following assumption (H) is satisfied 



(H): // d\x)d^{y) [ \m\\{77 



+ \xk + tk\)il + \yk - tk\) 



dt 



dxdy < oo. 



then 



rf+2Q+2/3 



(Qn-EiQn)) 



L?(0)L2(0) 



d{x)d{y) / g{t)H{x + t)H{y-t)dtdW{x)dW{y), (7) 



n^hereH{z)^UU'^- 

It seems that (H) is almost necessary for the convergence, since it guarantees 
that the integral in fT]) exists. 

The proof is relegated to Section 12.41 The main tool consists in rewriting 
Qn as a double stochastic integral. Then the convergence in quadratic mean is 
deduced from a simple convergen ce in L^(R^'^). Note that c ontr ary to the proofs 



for no n central limit theorems in iTerrin and TaqqiJ (Il99p|) andlDoukhan et al 



1996), we do not apply the scheme of convergence of iDobrushin and Maior 



1979I ). especially their Lemma 3 which involves a spectral measure convergence 



assumption. 



2.3 Some class of models satisfying the hypothesis (H) 

The s etting of random fields allows many kind of dependencies (see iLavancier 
1 20051 )). Apart from its intensity, the dependence can be isotropic or can occur 
all over several particular directions, depending on the form of the spectral 
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density. Since we did not want to restrict ourselves to one particular case, 
Theorem [2] involves a general integrability condition, i.e. assumption (H). We 
check in this section that this hypothesis is not too restrictive. Indeed, we prove 
that in all the preceding studies about quadratic forms of Gaussian fields, this 
condition is fulfilled. Moreover, it allows to extend to new cases, mainly when 
the dependence is not isotropic. 

As a consequence of the following Lemmas and according to Theorem [21 a 
non-central theorem for Q„ arises for all the cases treated in this section. 

Lemma 3. Assume that there exists some positive constants c and c' such that 



X, 



i5(^)i<c'nti 



x. 



\2l3/d 



(8) 



// moreover a > —d/2, jS > —d/2 and a + /3 < —d/A, then (H) is fulfilled. 



Remark 1. The result of our Theorem\M in the setting of Lemma\^ i s the same 
as the convergence stated in dimension 1 bu lTerrin and TaqqA ^199d ) (Th. 1). 

The assumptions in Lemma\^ allow the filter a{x) to he isotropic, i.e. equiv- 
alent at zero, up to a consta nt, to This is the hypothesis done in Theorem 
5.1 of lDoukhan et cH. 1 199d ) where, moreover, it is assumed (3 — 0. 

Proof. We have to check that the function defined on {M.'^)^ by 

~a\xra^{yMt)\\m\ 



(x,y,s,t) ^ 



nfc=i(l + + tk\){i + \yk - tk\){i + \xk + sfe|)(i + \yk - sk\) 



is integrable on (R'*)'*. From the hypothesis of Lemma O it is enough to prove 
that for all fc = 1 , . . . , d 

{tk,sk,xk,yk) ^ (1 + + tk\)il + \yk ~ tk\){l + \xk + sk\)il + \yk - sk\) 

is integrable on M^. Such integrals are studied in Lemma 1 of iTerrin and Taggu 
|1990) : Under the assumptions of Lemma[3]on a and /3, this integral is finite. □ 



When the functions a and g involved in Theorem [5] do not satisfy condi- 
tion (H) may be investigated thanks to power counting theorems (cf. Theorem 
2 in lTerrin and Taggu The following lemma focus on a particular sit- 

uation in dimension d = 2: When d admits one or two independent lines of 
singularities. 

Lemma 4. Assume that d ~ 2 and that g follows the same conditions as in 

m- 

If, for p^q, 

d{xi,X2) = \xi +pX2\°'''\xi+qX2\"''', 

and if ap > —1/2, > —1/2, (3 > —1 and ap + aq + /3 < —1/2, then (H) 
is fulfilled. 
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Proof. Checking (H) is equivalent to prove that 

f \Xl+pX2\^"-\xi+qX2\^"^\yi+py2\^''''\yi+qy2\^''^\tl\^\t2\^\si\^\S2\'' ^^^^^^^ 

(9) 

is finite. We apply Theorem 2 of iTerrin and Taaaul ()l990l) . Let us introduce 
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some notations. The integral above can be written : 
|Li (u) 1 I L2 (z.) I L3 H 1 IL4 (") P"' X 

\L5iu)f\Leiu)f\Lr{u)f\Lsiu)f Y[{1 + \Lkiu)\r'du, (10) 

k=9 

where u — (xi, X2,yi,y2iti,t2, Si, S2) and the Lfe's are the linear functionals 
involved in ([§]). For instance Li{u) = Xi +px2, L^{u) — ti, Lq{u) — X\ -\-t\. 

Let T — {Li, . . . , Lie} and let 7^ be the exponent associated to in PU)) . 
For instance, 71 = 2Q;p, 75 = /3, 79 = —1. Consider now the subsets W C T 
such that span{W) DT — W. A subset W is said padded if any G 14^ is a 
linear combination of the Li's in — {Lfc}. 

The integrability of ([9]) near i s obvious since 2q :p > — 1, 2aq > — 1 and 
P > —1. According to Theorem 2 in lTerrin and Taggu ( 199Cll ). the integrability 



at infinity is achieved if for every padded W considered above but T, doo{W) := 
rank{T) - rank{W) + J2t-w Tfe < 0- 

The maximum value for doo{W) is obtained with W — {Lg, . . . , Lie}. In this 
case doo{W) = 8 - 6 + 4ap + 40;, + i/S. This leads to + a, + /3 < -1/2. □ 

2.4 Proof of Theorem M 

Pr oo f. Le t Z be a Gaussian random spectral measure associated to the measure 



H. iMaioij (|1981D defines 

f{x,y)dZix)dZ{y) 



for all / G where denotes the space of functions /: M'' x E'' — >■ C such 
that f{—x,—y) — f{x,y) and J \f{x,y)\'^dp{x)dp{y) < 00. The second order 
moment of this integral satisfies, for all / in iJ^, 

E (^JJ fix,y)dZ{x)dZ{y)^ < 2 J^Jf{x,y)fdti{x)dfi{y). (11) 
And the so-called Ito formula follows: For all /i and /2 in iJ^, 



h{x)f2{y)dZ{x)dZ{y) = I h{x)dZ{x) I f2{y)dZ{y)- / h{x)f2{x)d^i{x). 

(12) 
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Now, let us rewrite Qn as a double stochastic integral. According to ([T]), 
Qn ~ E{Q„) ^^Y.J2 9k~i {XkXi - r{l - k)) 



k&A„ leA„ 



fc6A„ IgA„ 

From the Ito formula 
XkXi - r{l - fc) 

a{x)e'<^''-'>dW{x) I a{y)e'<^^y>dW{y) 



dt. (13) 



£2 



a(a;)a(y)e 



i(<fe,a;> + <;,jf>) 



i</-fc,a;>„2 



a {x)dfi{x) 



dW{x)dW{y), 



where /x, the spectral measure of e, is proportional to the Lebesgue measure A. 
Assume that g{—x) — g{x) and a{—x) = a{x), then the function 



{x,y) ^ a{x)a{y) 



\ I git) V e^<'='-+*> V e'<''y-'>dt 



belongs to Ha- Therefore we can rewrite (fT3|) as 



(14) 



a{x)a{y) 



£2 



-4 / g{t)Hn{x + t)Hn{y-t)dt 



(iM^(a;)(iT/F(?/), 
(15) 

where i?„(i) = S,g^^e»<'=^*>. 

In (|15|l . we make the change of variables x — > x/n, y y/n, t — > t/n. 
Since the Gaussian measure W satisfies, for all Borelian set A, W{n^^A) = 
n~^/'^W{A) and since a and g are homogeneous, we get 

n''+'"+'^(g„-i?(Q„)) = 

a(n)(a;)a(„)(2/)ii Li '4)n{x,y)dW(x)dW{y), 



where 



^ + t\ 1 



dt. 



As a consequence, according to ([TT|). it suffices, for proving ([7]), to show that 
the following integral tends to zero as n — >■ oo: 



a(n)(a;)a(„)(?/)Li (^^^ Li (^^^ 5'„(a;, y)-a(„)(x)a(„) (2/)Li(0)^'(a;, y) 



dxdy 
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^{x, y) - L2 (0) / mH{x + t)H{y - t)dt. 



where 



From the decomposition A^B - CD = - C)B + {B - D)C, the i^-norm 
above is lower than the sum 2(/i + I2) where 



h 



and 

I2 = 



a(n){x)a(ri){y)'^n{x,y) - a(n){x)a(n){y)-^{x,y) 



dxdy 



L 



n 2 



a(„)(a;)a(„)(y)*(a;,?/) 



dxdy. 



The fohowing lemma will be useful. 
Lemma 5. 

(j) V z = (zi, . . . , Zrf) e ni?, 



lim 

n— >oo 



(ill) For a.e. z e R'' 

Proof of Lemma\^ Since for all j, < mr 

=n 



j=i V I 



1 e"^- - 1 
-e'-— 

n e*^ - 1 



n 



sin(zj/2) 



z,/2 



Zj/2n 



sin(zj/2n) 



Similarly, for all z e R'', 



i^wi=n 



e»^. - 1 



n 



sin(zj/2) 



z,/2 



Finally, for proving (Hi), suppose first that d — 1. If z 7^ 0, 



n \n/ 



1 , . e" - 1 e'^ - 1 

-e*" — 

n e'" — 1 iz 



2 


e "iz 








n(e'^ - 1) 
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The norm in the right hand side term is equivalent to \z\ln when n goes to 
infinity, hence it tends to 0. In dimension d, (iii) is proved by induction thanks 
to the decomposition AB-CD = {A- C)B + (S - D)C. □ 

Let us first prove that Ii asymptoticaUy vanishes. Since Li is bounded, it 
suffices to prove the convergence to of 



hi 



<i>i")(i,x,y)di 



71 E 



dxdy, 



(16) 



with 



and 



/l2 = 



aL(a:)aL(y)L2(0)2 



nE 



dxdy, (17) 



with 



—.Hn f ^) \Hn ] -H{x + t)H{y ~ t) 



In both III and /12, the 27r-periodicity of g, iJ„ and H allows us to reduce 
the domain of integration nE (with respect to t) to 

nDx^y = {\x -t\ < nn} n {\y + t\ < mi} n nE. 

Therefore, (i) of Lemma [5] can be applied and since L2 is bounded, 

d 



^i^^\t,x,y)dt 



nE 



< C 



l5(n)(i)in 1 



fe=l 



A 



1 



\xk + tk\ 



1 A 



1 



IVk - tk\ 



dt, 



where c is a positive constant. Hence, 



hi < c' 



c' ff -a'ixya'iy) f |g(t)| J] ^ 

JJu^" [Jr" j^J-^{l + \xk + tk\){l + \yk-tk\) 

Besides, according to (i) of Lemma [SJ 



dt 



dxdy. 



Ke-{x, j/)II„d..„ (ti)lI„D..„ {t2)dln) i^)aln) {y)^i''\ti , X , 2/)$l"Ht2, X, y) 



<c~a\x)d\y)\g{ti)\\~g{t2)\ 



i2(0) 



^2(0) 



From the continuity of L2 at 0, this term tends to zero for any fixed (x, y,ti,t2) S 
M^"^. Therefore, thanks to assumption (H), the Lebesgue's dominated conver- 
gence theorem applies and lim„_>oo hi = 0. 
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The convergence of I12 is proved similarly. From (i) and (ii) of Lemma [SJ 
we have 



nE 



<c/ \~9in)it)\l[(lA . i.i VlA ^.] dt 

JnD,_y fJl\ \xk+tk\J \ \yk-tk\J 



and 

h2<c [[ ~a\xya^{y) 



k=l 



{l + \xk+tk\){l + \yk-tk\) 



dt 



dxdy, 



where c and c' are positive constants. 

Besides, for almost every {x,y,ti,t2) £ M^'*, according to (iii) of Lemma [SJ 
Ke^ (x, y)JLnD^,yiti)TLnD^,, (i2)5(„) (a;)a(„) {y)<^>l^^^ {ti,x, y)$l^^) {t2,x, y) ^^-^ 0. 

The Lebesgue's dominated convergence theorem applies thanks to (H) and /12 
tends to zero when n ^ 00. 

It is easy to see that I2 tends similarly to zero. □ 



3 Convergence of the empirical auto-covariance 
function 

We present an application of the preceding theorems to the asymptotic law of 
the empirical covariance function in U^. Indeed, for h e Z"*, 



i£A„ 



is a particular case of ([3]) with g{t) = {2tt 



,—d'i<h,t> 



We will consider further 



where Xn = n^'^ J2ieA -^i ^^"^ denote in the following r{h) ~ E{XiXi+h)- 
We are dealing with linear fields as ([!]). In dimension 1, this framework 
provides a dichotomy in the asymptotic behavior of f (fe) depending on whether 
the spectral density of X belongs to or not (see iHoskina ( 19961) ). We prove 
later on the same kind of division in dimension d, although some kind of in- 
termediate behavior can arise. Let us note that thi s dichotomy does no longer 
hold in d = 1 when X is not linear, see for insta nce Giraitis et al. ( 200Cl( ) in the 
case of a LARCH process or Giraitis and Taggu ( 1999( ) for the consequences on 
the Whittle estimator. 

When / e L^, we prove that f{h) and f{h) follow the same central limit 
theorem. This is the object of Proposition [6l 
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When / ^ L^, the asymptotic behavior of f{h) and f{h) may differ and we 
focus on f{h) which is a proper quadratic form. In this setting, the asymptotic 
law comes from Theorem [5] where /3 = 0, provided condition (H) is satisfied. 
Therefore, the normahzation and the hmit in law depend on the filter a. We 
summarize in Proposition[7]and in the beginning of section the two situations 
already studied in Lemma |3] and Lemma 0] before. 

But Theorem [5] does not apply in the example of Lemma 2] when UpUq = 
and /3 = 0. This corresponds to the particular situation when the long memory 
occurs only along one direction in dimension 2. In this case, we study the 
asymptotic behavior oif{h) in PropositionlHl It appears that a non-central limit 
theorem holds in the sense that the normalization is not n'^/'^. Yet, contrary 
to the classical non central limit results for f{h) (see references therein), the 
limiting law is Gaussian. This shows that the covariance series of a linear field 
may be asymptotically gaussian even if its spectral density does not belong to 

L2. 



3.1 General results in dimension d 

Let us first present the central limit theorem for f{h) and f{h) when f ^ . 

Proposition 6. Let X be the linear field defined by (|7J). 

If f e i^ then, for all h £ Z'^ , n'^/^{f{h) - r{h)) and n'^/^{r{h) - r{h)) 
converge both in law to a Gaussian random variable with zero mean and variance 
2(27r)''/|^, where /|^ stands for the {2h)-th Fourier coefficients of f^ as defined 
in 

Proof. The central limit theorem for r{h) follows from Theorem [T] where p = 2, 



-oo and g(t) = (27r) 



-d^i<h.t> 



For r(/i), we prove that almost surely, n'^^^{f{h) — r{h)) — o(l). Indeed, 

1 



,3d/2 



Using representation ([T]) of X, this term is equal to 



1 



,3<i/2 



i<h,x> 



E 



dx 



1 



,d/2 



f{x)e<'''^>\{F^{xj) dx 



where Fn denotes the Fejer kernel on [— 7r,7r]. Now, from the Cauchy-Schwarz 
inequality, this last term is lower than 



1 



^d/2 



\ 



P{x)X{Fn{Xj) dx. 



Since p G L^, the Lebesgue's Theorem implies that the integral above is a 
o{nf'). Therefore, n'^^'^{f{h) — f{h)) = o(l) and the central limit theorem for 
f{h) is inherited from the one for f{h). □ 
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Let us recall that the notation a(i) is defined at the beginning of Section 
12.21 The following proposition provides central and non central limit for f{h), 
depending on the parameter a, when condition ([T5)) below is fulfilled. This 
framework includes isotropic and anisotropic models. The result shows the same 
kind of dichotomy than in dimension d = 1. In section 4.2, some anisotropic 
models that do not follow this dichotomy are presented. 

Proposition 7. Let X be the linear field defined by |IJ). 

Assume that for all x ^ E a{—x) — a{x) and that a{x) = a(i)(a;)Li(x), 
where a is a homogeneous function of degree a > —d/2 and Li is a bounded 
function, continuous at zero and non-null at 0. If 

d 

\a{x)\<c\[\x^\''''^ (18) 

i=l 

for some positive constant c, then for all h G Z'', 

• if a > -d/A, 

n'/'ifih) - r{h)) A (O, 2i2nrfi,) , (19) 

where is the {2h)-th Fourier coefficients of f^ as defined in 

• if a < -d/A, 

n'^+^"{f{h)-r{h)) Ll{0) [[ &{x)a{y)H{x + y)dW{x)dW{y), (20) 
where Hlz) = jt- , 

Proof. In the case a > —d/A, the convergence result is a consequence of Propo- 
sition [51 

For a < —d/A, Lemma |3] applies since /3 = and condition (H) in Theorem 
[2] is fulfilled. Let us justify the simplification of the limit in (pO)) . Here 

g{t) = (27r)-'^e^<'''*>, 

so g{t) = 1 and ^2(0) = {2tt)~'^. The simplification comes from the main term 
J12 in the proof of theorem [51 where we use the identity 

JnE " \ n J \ n J n'^ \ n J 

The pointwise convergence of this last term relies on {iii) of Lemma [SI and an 
application of the Lebesgue's theorem concludes the proof of (|20p. □ 
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3.2 Some anisotropic examples in dimension d = 2 

Starting with the anisotropic case studied in Le mma [H we confirm th a t new 
limiting results can be obtained as suggested by iDobrushin and Major ( 1979f ) 
(Remark 4.2). If we suppose that a{x) = d(^i){x)L[x), where L is a bounded 
function, continuous and non-null at and 

d{xi,X2) = \xi +PX2\"''\X1 +qX2\°''', 

with ap > — i, aq > — ^ and (p, q) G R^, then, under the assumptions of Lemma 

where Z is defined as the limit in (I^U)) . 

This convergence is a simple application of Theorem [21 Condition (H) is 
fulfilled thanks to Lemma |4] and the simplification of the limit holds for the 
same reasons as for Proposition [71 

The assumptions of Lcmma[4l in the case when /3 = 0, imply the existence of 
two lines where the spectral density is unbounded. When ap — or aq — (or 
p = q), that is when the long memory occurs along only one direction, Lemma 
[4ldoes not imply (H) when /3 — 0. This case provides a new limiting behavior 
for f{h) as stated in Proposition [51 

Before stating this result, we deduce in the following lemma the covariance 
structure of X when the spectral density is unbounded along a line that goes 
through the origin. We say that the process has a long memory along one 
direction if its covariance function is not summable along this direction. 

Lemma 8. Let p e M and for all {xi,X2) E M^, 

f{xi,X2) = ^/(^^i +P^2), 

where f is an even, non-negative, and 2TT-periodic function on R. 
Let us denote, for all (hi, 112) 6 Z^, 

a{hi,h2) = I e*(''i"i+''^"^V(a:i,X2)dxidx2 

J[-7r,7r]2 

and 

/TT 
{x)dx. 
-TT 

Then, for all (/ii,/i2) £ Z^, we have 

, s i^{hi) ifh2 ^phi, 

j,Hiffiii_glM^(/,^) otherwtse. 

The proof of this lemma is given at the end of this section. 
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Proposition 9. Let X he a stationary Gaussian process in dimension d = 2. 
Let us suppose that its spectral density is 

where p £ Z and where f, defined on [— 7r,7r], is a spectral density in dimension 
d = I. Assume moreover that for — 1/2 < a < and for all x G [— 7r,7r], 
f{x) = where L is a bounded function, continuous at zero and non- 

null at 0. Then, 

• ifa> -1/4, ^ 

n{f{h) - r{h)) A U (O, 2(27r) VI,) , (21) 

where /|^ is the {2h)-th Fourier coefficients of f^ as defined in 

• if a< and if h2 ^ phi, 

n(4"+3)/2(f (/j) _ r{h)) A A/- (0, alp) , (22) 
where <j'^ p — limn^oon'^"^^Vo,T{f{h)). 

Proof. In the case a > —1/4, the resuU is a consequence of Proposition [B] Let 
us focus on a < —1/4. We restrict the proof to the case p > since p < can 
be treated in the same way. 

We p rove the result thanks to a central Umit theorem for triangular arrays 
stated in Romano and Wolf (I2OOOI) . 

An alternative proof could be to start from the representation of f{h) in 
terms of a double stochastic integral as in (|15p , then to use the n ecessary and suf- 
ficient cond itions for the convergence of s uch an integral given in lNualart and Peccati 
12005.) and Nualart and Ortiz-Latorre ( 20081 ). But checking these conditions 
leads to hard computations, which would not simplify or shorten the present 
proof. 

So, let us take advantage of the natural triangular array structure that arises 
in our setting. Indeed, let Yi-^^i^ — Xi-^^pi^^^i^ and let i = ii, j — 12 — pii, we 
have 



^ n n ^ n—p 

f (/ll, /I2) = ^ ^ ^ -'^Il,l2-'^il + ?ll,i2+/l2 = ^ ^ ^ Kj^jFi + Jl^ 

ii=li2 = l j =—pn+lieBj 

where Bj = {i \ 1 A ^—^ < i < V n}. Therefore r{hi, ^12) is the triangular 



array 



where 



{p+l)n—p 
f{hi,h2) = ^ YnJ 

Yn^j — ^ ^ ^ Yi^j—pnYi-\-}ii,j—pn—phi+h2- 
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Let us summarize the properties of Ynj- 
From Lemma [51 we have 

^, , ^ }r{hi) if h2 =phi, 
r{hi,h2) = < . (23) 

I otherwise. 

where f is the covariance function associated with /. Consequently (Yij) is a 
zero mean Gaussian process such that 



P{hi) if/l2=0, 

otherwise. 



The process F„j , viewed as a function of ((Fi j_p„)igz, {Yi+h,j-pn~phi+h2)iei.)^ 
is thus a (/i2 — p/ii)-dependent process. Moreover, we can compute the moments 
of Ynj thanks to the representation of the moments of Gaussian variables in 
terms of Wick's product. Since we have assumed h2 ^ phi, we obtain 

E{Yr,^j)^0, 



E{Yn.j-i Y,i j2 

and 



otherwise 



6n ^ ^ r(i2 - ii)r(i3 - i2)f(i4 - i3)f(i4 - ii). 



We are now in position to apply Theorem 2.1 in lRomano and Wo la ()2000D which 
gives sufficient condition for the convergence in law of a triangular array of 
m-dependent random variables to a normal distribution. Following the same 
notations as in this theorem, we choose 6 — 2 and 7 = and we look for A„, 
Kn and L„ such that 

. E{Yl^)<A^ for all J, 

• Var [J2'i=a^^ ^n.j) < k Kn for all a and for all k > (/12 — phi), 



• Var (Eg'^""' Yn,j) > {{P + l)n ~ p) 



According to iRomano and Wolj (l2000l) . the convergence in law holds whenever 



KjL.n = 0(1) and A„,/Lg = 0(1). 

From Theorem 2.24 of IZvgmundl ( 1959f ). we have r{h) Cah~^"~^ when 



/i — 00 where is a constant depending on a and L{0). 
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?.a-8\ 



Since {i € Bj) ^ {1 < i < n), an integral test leads to A„ — 0{n~ 
and Kn = 0(n~'*"~^). For L„, let us compute directly the equivalent of the 
variance. From the decomposition in terms of Wick's product, we obtain 



(j9+l)n-p 



n 



,-4 



n ■ ^ r{i2-ii,j2-ji)+r{i2+hi~ii,j2+h2-ji)r{i2~ii~hi,j2~ji-h2)- 

In view of (|23p and since /12 7^ phi- most of the terms in the sum above 
vanish. We finally obtain 

((p+l)n-p \ 
i=i / (n,i2,ii)ec 

where C = {1 < ii, 12, ji ^ n, 1 < ji +p{i2 — *i) < JT-}- Let fc = 12 — ii, we have 
C = {l<j<ri, |fc|<n. — 1, 1 — pk < j < n — pk} 

_ 1 

<A:<— 1, 1 — pk < j < n 



|o < fc < 


n-l 






. P . 



So, 



(p+l)ri— p 

i=i 





71-1 


H- 













■{0) + 2n-^ E r^{k){n~pk). 



k=l 



The latest sum involves positive terms which are equivalent to c^k^'^°'^'^{n~pk) 
when fc — !• 00. As a consequence, Var (^X^j-S^"*" ^^n,i) ~ "'a p"-~^"~^i where 
tr^ p is a positive constant. This lead s in particular to L „ = 0{n~'^°'~'^). 

The conditions in Theorem 2.1 in Romano and Wolf ( 2000l) are fulfilled and 
the convergence in law holds. □ 

Proof of Lemma\^ When p = 0, the result is obvious. Let us assume, without 
loss of generality, that p > 1. 



a{hi,h2) 



1 

2^ 



g*/«i(:.i+px2) + px2)e''^'''-P''''>''''dxidx2 



Let the change of variables u — xi + px2 and v — X2 ■ 



a{hi,h2) = — 



1 Ap+i)^ f r^^^ \ 



-{P+1)7T 
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Let us first suppose that h2 ^ phi. When p > 2. the above domain of 
integration can be cut up as follows (the case 1 < p < 2 is not detailed but can 
be treated similarly): 

where 

-(p-1)t 



h = -i [ e'^i''/(u) U^^-'-P''^^'^ - e-»(^2-p/»i)7rA 

/•(p+l)7r 

J(p-i)v ^ ' 



'(p-i) 

and 



J-(p-l)7r ^ ^ 



Some trigonometric computations lead to 

f (p+l)7r 

'(p-l)7i 



7i = 2 / /(u) I sin(/iiM + (/i2 — phi)^) — sin (/12 — phi) — 

J(p-l)w \ \ P P 



du. 



h = 4(-l)"'iszn / /(u) cos du. 



P J Jo \ V 

Now, let s = M — \jp\ and e= p — [pj , 

/(s + [pJtt) (sin(/iis + /i2 LpJtt + (/12 - p/ii)7r) 

-TT+eTT 



. / h28 , LpJtt ,, , ,7r\\ , 

■sm h ft2 (/i2 - phi)— ds, 

P P Pjj 



1, = 4(-I)^. Sin [^) (J^"'-'" /(„) c™ 

+ /(. + LpJ.) cos (^ill + ft.Mil) rf.) . 

The domain of integration in Ii can be split into — tt < s < tt and s G 
TT, — TT + eTr] U [tt, TT + rtt] . From the 27r-pcriodicity of /, this is easy to check 
that, when summing up Ii and I2, all the integrals involving e in their range of 
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integration sum up to zero. Hence Ji + I2 reduces to 



/(s + [p\n) (sin(/iis + /i2bj7r + (/i2 - phi)Tr) 

. ( h2S , [pJtt , n 1 7 

—sm h /i2 (n.2 — ph'i)— as 

\ P P P, 



/(u) cos — ) dw. (24) 



P 



The latest integral above is 

LpJ-2 

/ .f(u+ jTT)cOS i h2 

P 



LpJ-2 ^ ^ 

J=0 ^0 



)co.s ( /^s ^"^^""^ ]du (25) 



and it is handled according to the parity of [p] and j. Since / is a 27r-periodic 
function, f{u + jtt) = f{u) when j is even and f{u + ji:) = f{u — vr) when j is 
odd. 

When Ip\ is even, the sum above is then 



/ /(") 2^ cos U2 



du 



P 



h2u\ sin (^h2^- 112 



sin h2^-h2^ 



f{u) COS ( ) J r — —du 



When plugging in this latest result in ([M)) , it simplifies and Ii + 12 becomes 

r ~ 

f{s)sm{his + {h2 — phi)Tr)ds ^ 2sin{{h2 — phi)TT) / f{s)cos{his)ds. 



This proves the result of the lemma for /12 ^ phi in the case [p\ even. 
When [pj is odd, the sum p5[) is 



LpJ-3 

^27r 



I fiu)±C0.(^h2^-^^)du^ 



J=0 

f.27r 



COS ^ 






p 2 J 


sin ( 









/(u) sm ( — ) ^ j^—^ '-du. (26) 

Now, in we split the domain of the first integral into — tt < s < where 
f{s + [pjvr) = f{s + tt) and < s < tt where /(s + [pJtt) = /(s — tt). We apply 
respectively the change of variables s ~ s + t: and s = s — vr. This allows to 
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exhibit the integral 2 J^^ f{s) sm{his + (/12 — phi)-K)ds. With the help of (j^ . 
some trigonometric computations show that the remaining terms coming from 
this change of variables simplify with the remaining term in (j24p . 

Therefore, when /12 ^ phi, the result of the lemma is proved for all p. 

The proof when /12 — phi is simpler and it can be conducted in the same 
way. □ 
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